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o ' 

, We discuss the transition from quasi-circular inspiral to plunge of a system of two non- 

■ rotating black holes of masses m± and ni2 in the extreme mass ratio limit m\m,2 <C 

O ' (mi + n%2) 2 . In this limit, we compare the merger waveforms obtained by two differ- 

ent methods: a numerical (Regge-Wheeler-Zerilli) one, and an analytical (Effective One 
Body) one. This is viewed as a contribution to the matching between analytical and 
numerical methods. 
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1. Introduction 



I/"") , The last months have witnessed a decisive advance in Numerical Relativity, with 

different groups being able to simulate the merger of two black holes of compara- 
ble masses. 1 Since such binary black holes systems (of a total mass^ 30M Q ) are 
\Q \ believed to be among the most promising sources of gravitational waves for the 

ground based detectors like LIGO and VIRGO, this breakthrough raises the hope 
O , to have, for the first time, a reliable estimate of the complete waveform by join- 

ing together Post-Newtonian (PN) and Numerical Relativity results. We recall that 



(N 



■ PN techniques have provided us with high-order results for describing the motion 



>< 



and radiation 2 of binary systems, and that further techniques have been proposed 
for resumming the PN results, 3,4 thereby allowing an analytical description of the 
gravitational waveform emitted during the transition between inspiral and plunge, 



and even during the subsequent merger and ringdown phases. We now face the 
important task of constructing accurate complete waveforms by matching together 
the information contained in Post- Newtonian and Numerical Relativity results. We 
view the present work as a contribution towards this goal (see also Refs. 5 ' 6 ). 

The present work belongs to a scientific lineage which was started by Regge and 
Wheeler, 7 Zerilli, 8 Davis, Ruffini, Press and Price 9 and Davis, Ruffini and Tiomno. 10 
Refs. 9 ' 10 studied the gravitational wave emission due to the radial plunge (from in- 
finity) of a particle into a Schwarzschild black hole, as a model for the head-on 
collision of two black holes in the extreme mass ratio limit. Here we consider, for 
the first time, the transition from the quasi-circular adiabatic inspiral phase to the 
plunge phase in extreme-mass-ratio binary black hole systems. We shall be able 
to do that by getting round a limitation of the original Regge-Wheeler-Zerilli test- 
particle approach: that of requiring the test particle to follow an exact geodesic 
of the Schwarzschild background. We bypass this stumbling block by using an im- 
proved form of PN theory: the Effective One Body (EOB) approach to the general 
relativistic two body dynamics. This approach has been recently proposed to study 
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Fig. 1. Left panel: Gravitational waveforms (real and imaginary part) for I = m = 2 generated 
by a plunge from an initial separation r = 1M. Right panel: Comparison between the corre- 
sponding gravitational-wave phase and that obtained from analytically matching a 3PN improved 
quadrupole-type formula to a superposition of quasi-normal modes. 



the transition from inspiral to plunge in the comparable-mass case. It describes the 
dynamics of a binary system in terms of two separate ingredients: (i) a Hamilto- 
nian H^ob{M^ A*) describing the conservative part of the relative dynamics, and (ii) 
a non-Hamiltonian supplementary force .Feob {M, ^i) approximately describing the 
reaction to the loss of energy and angular momentum along quasi-circular orbits. 
[Here, M = mi + m,2, and /i = m,\milM\ The badly convergent PN- Taylor series 
giving the angular momentum flux is resummed by means of Pade approximants. 3 



2. Results 

We summarize here our tools and main results (see Refs. 12,13 for more informa- 
tion). We describe, in the extreme mass ratio limit, the relative dynamics of the 
binary system of two non-rotating black holes by that of a "particle" (of mass 
[i <C M) that follows a non-geodesic inspiral driven by radiation-reaction (until a 
quasi-geodesic plunge) on a (quasi-) Schwarzschild background of mass M. We then 
use the multipolar Regge- Wheeler- Zerilli perturbation theory around such a black 
hole to compute the gravitational wave emission in the approximation /i <C M. 
Einstein equations then lead to two (even/odd) decoupled wave- like equations (with 
source terms 5^°' oc /i) for two (even/odd) master functions - 11 

Fig. 1 shows the I = m = 2 waveform generated by a binary system with initial 
relative separation r = 7M 0. After a "chirp-like" increase of the frequency and 
modulus during the inspiral, the latter reaches a maximum when the particle crosses 
the light ring (r = 3Af ) and eventually the waveform decays in a quasi-normal mode 
ringing phase. 



a Note that the other multipoles give a non-neglible contribution, due to the strong asymmetry of 
the system. 
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The so-constructed "exact" numerical waveforms are then compared with semi- 
analytical ones, constructed within the EOB framework and philosophy. We recall 
that the basic idea of the EOB framework is to produce quasi-analytical waveforms 
by patching together a quadrupolc-typc waveform during the inspiral and plunge 
to a QNM-type waveform after merger. We give an example of this numerical- 
analytical comparison in the right panel of Fig. [TJ The gravitational wave phases 
(fy^ obtained by two different methods, one numerical and the other semi-analytical, 
are compared. In the first method, 4^ is given from time-integration of the instan- 
taneous gravitational wave frequency obtained as Muj 2 ^ = -S>(*^/*22 )■ Tnc 
second method computes an approximate waveform in the following way: before 
crossing the light ring one uses a (Pade resummed) 3PN-improved quadrupole-type 
formula to compute the waveform from the EOB dynamics. After crossing the light 
ring, the previous quadrupole-type signal, taken in the quasi-circular (QC) approx- 
imation, is matched to a superposition of the first five QNMs of the black hole. 
Then one computes the phase of the matched analytical waveform by integrating 
the corresponding instantaneous gravitational wave frequency. The maximum dif- 
ference between the "exact" phase and the "effective one body" phase turns out to 
be less than 3% of a cycle (see inset). 
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